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The exact solution of the boson pairing hamiltonian given by Richardson in the sixties is used 
to study the phenomena of level crossings and quantum phase transitions in the integrable regions 
of the sd and sdg interacting boson models. 



One of the most fruitful features of the Interacting Boson Model (IBM) [1] of nuclei is the existence of three 
dynamical symmetry limits. Each represents a well defined nuclear phase, providing analytically the exact eigenstates 
of the system and offering a unique tool to deeply understand the physics involved. 

A quantum system has a dynamical symmetry (DS) if the hamiltonian can be expressed as a function of the Casimir 
operators of a subgroup chain. A direct consequence of this definition is that a system exhibiting a DS is quantum 
C\| \ integrable and analytically solvable. The concept of quantum integrability (QI), though sometimes associated with 
,i— > ■ that of dynamical symmetry, is more general, however. It can be stated as follows: A quantum system is integrable 
if there exists a complete set of hermitian, independent and commuting operators (constants of motion). Clearly, if a 
\^ , quantum system has a DS, the Casimir operators of its subgroups play the role of the constants of motion, fulfilling 
' the definition of QI. But, by no means must a system have a DS to be quantum integrable. 

Perhaps the most striking manifestation of quantum integrable systems is the absence of level repulsion in their 
t-H ■ energy spectra. This can be clearly seen in the IBM by looking at the transitions between the three dynamical 
symmetry limits. Using the Consistent-Q hamiltonian [2], 

H = xn d + ^-Q x ■ Q x (1) 
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the transitions between the dynamical symmetries 



In Fig. 1A, we plot the + energy levels associated with the transition from SU(3) (x = 0, x = — \/7/2) to 0(6) 
(x = 0, x = 0) and with the transition from 0(6) to SU(3) (x = 0, \ — y/7/2). The results are shown as a function 
of x and for N = 10 bosons. The level repulsion between pairs of energy levels is present everywhere, except at the 
X — point where the 0(6) DS is realized. At this point, due to the absence of level repulsion, level crossings are 
£h ' allowed. In the inset we amplify the first of the two allowed level crossings. 

In Fig. IB we show the analogous results for the + states of a system of TV = 10 bosons along the transition from 
SU(3) (x = 0, x — ~ v7/2) to U(5) (x — 1), as a function of the parameter x. The level repulsion phenomenon 
appears any time two levels are close enough. The inset amplifies one of the closest two level approaches, showing 
that level repulsion prevents the crossing. The level repulsion between the ground state and the first excited + state 
leads in the thermodynamic limit to a non-analytic behavior of the ground state energy, defining a critical point 
characterized by a first-order phase transition [3]. 

An interesting question is what is the structure and origin of the phase transition from SU(3) to SU(3). As noted 
earlier, this transition proceeds through the 0(6) DS. As also noted earlier, level crossings are permitted at that 
critical symmetry point. Thus, if the SU(3) to SU(3) transition is indeed first-order, it is interesting to ask how this 
can come about in the absence of level repulsion at the critical point. Further study on this issue is clearly required. 
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The transition from 0(6) to £7(5) is described by the hamiltonian (1) by setting x = and varying x from to 1. 
Earlier investigation of the properties of the spectrum for this leg of the Casten triangle made clear that the system 
is integrable everywhere on the path, even though there is no global DS [4-6]. To facilitate discussion of this region 
of parameter space, it is useful to digress for a moment and to discuss a class of integrable boson pairing models. As 
we will see, the 0(6) to £7(5) leg of the Casten triangle falls within this class of models. 

The models of interest are among three classes of integrable models for fcrmion and boson systems [7]. Though 
integrable, these models do not in general relate to a DS. The models we will be discussing are based on the SU(l, 1) 
pair algebra 



(2) 



where the operators b[ m (£>/ TO ) create (destroy) a boson in state Im. In terms of the generators (2), the constants of 
motion of the rational class are 



Rt = K? + 2g ]T 
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It can be readily confirmed that the operators Ri are hermitian, are independent (none of them can be expressed 
as a function of the others), and mutually commute. Moreover, if the system has M single boson states I, there are as 
many operators (3) as quantum degrees of freedom, constituting a complete set of constants of motion. The pairing 
strength g and the set of M real numbers vi are free parameters. 

If a system is quantum integrable there should exist a unique basis of common eigenstates of the M operators i?;. 
It has been shown [7] that this complete set of eigenvectors can be formally written as a product of boson pairs acting 
on a subspace of unpaired boson states \v) = \v\, V2, • • • , vm) 



i*>=n e 

a=l \Z=1 




(4) 



where the n parameters e a that apply to each eigenstate are particular solutions of the set of n coupled nonlinear 
equations 



' 2vi-e a 9 e a -e p ~ 



(5) 
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The total number of bosons is N = 2n + J2i v l an d the eigenvalues r; of the R\ can be found in Ref. [7]. 
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Any hamiltonian constructed as a linear combination of the constants of motion (3), viz: H = 2^ ; eji?j (g,rj), 
commutes with them, and thus is diagonal in the basis of cigenstates (4). The hamiltonian eigenvalues are the same 
linear combination 2^ ; e/n (g,rj). In particular, the boson pairing hamiltonian that was solved by Richardson [8] 
is obtained by choosing the coefficients of the linear combination equal to the parameters inside the Ri operators, 
ei = rji. We will use the following parameterization of the pairing hamiltonian: 

H P = lH lb \mbi m + 1 -^- £ b\ m b\-b vm ,b Vm , (6) 

Ira ImVm' 

for which the eigenvalues are 

M 

e p = Y1 + £ e " ( 7 ) 
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With the limitation to s and d bosons, the hamiltonian (6) describes a transition from a spherical vibrational U(5) 
DS to the 7 unstable 0(6) DS. In Fig. 2 we show the corresponding + states for N — 10 bosons as a function of x. In 
the limit of x = 0, we are in the exact 0(6) DS limit; for x = 1 we are in the U(5) DS limit. As previously discussed, 
the behavior of the levels as a function of x does not show any sign of level repulsion. Since the 0(5) Casimir operator 
commutes with the hamiltonian (6), the cigenstates can be labelled by the boson seniority quantum number r. 

In the left panel, we plot all the + levels, with the seniority quantum number r specified on the right vertical 
axis. There are several level crossings in the figure, but they all correspond to pairs of levels with different seniority 
quantum numbers. The + levels with r = are displayed in the right panel. These levels evolve independently with 
the parameter x. The fact that the complete transitional region from 0(6) to U(5) is quantum integrable has been 
been previously noted [6]. From the six quantum degrees of freedom of the C/(6) dynamical group, four of them are 
taken into account by the Casimir operators of the 0(5) subgroup chain, the fifth is the conserved number of bosons 
N and the sixth is any linear combination of the Casimir operators of 0(6) and U(5), e.g. the hamiltonian, which by 
construction commutes with all the other constants of motion. We arrive to the same conclusion from the constants of 
motion given in (3). In an sd space there are two constants of motion, R s and Rd- The sum gives the boson number 
N and any other combination defines a hamiltonian interpolating between 0(6) and U{b). The states with unpaired 
bosons are completely classified by the 0(5) subgroup. 
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FIG. 2. + energy levels for a system of iV = 10 bosons as a function of the parameter x for the 0(3)-C/(5) transition. The 
left panel shows all + levels; the right panel only the r = levels. 



A non-analytic point in the thermodynamic limit due to level repulsion is precluded in the transition from 0(6) to 
?7(5). The only source of non-analyticity within an integrable region is a level crossing. However, as we see in Fig. 
2, states with r = do not cross. The exact solvability of the model implies that the energy is an analytic function 
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of the parameter x and only phase transitions of order greater than one are permitted. We may wonder, therefore, 
under what conditions are crossings between t = states possible. It can be shown that crossings between states 
with the same set of quantum numbers can take place when there are at least two parameter dependent constants of 
motion. The simplest example is a boson model with sdg bosons. In this case, there are three independent constants 
of motion, R s , Rd and R g , but factorizing out the boson number N as the sum of the three, we are left with two 
that are parameter dependent. An analysis based on the structure of the U(15) dynamical group of the sdg IBM 
leads to the same conclusion. The pairing hamiltonian (6) is a linear combination of the Casimir operators of two 
subgroups, those associated with ?7(14) and 0(15). Both subgroups constitute a DS of the 17(15) dynamical group 
and both have the group 0(14) as a common subgroup. In Fig. 3 we show the + , t = states for N = 10 bosons 
as a function of the parameter x for this model. As expected, there are no signs of level repulsion, but there are level 
crossings. However, there are no crossings with the ground state, which would have been evidence for a first-order 
phase transition, even though it is a finite system. 
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FIG. 3. sdg + ,r = energy levels for a system of N = 10 bosons as a function of the parameter x for the 0(15)-Z7(14) 
transition. 

In previous works [9,10], we showed that the rational class of integral models with repulsive pairing exhibits a 
quantum phase transition (QPT) to a symmetry broken phase with a macroscopic occupation of the two lowest 
single-boson states. This led us to conclude that repulsive pairing between bosons is a new and robust mechanism for 
enhancing sd dominance [10] in interacting boson models of nuclei. 

Within the Landau theory, a second-order phase transition is related to a continuous change of the order parameter 
from in the disordered phase to a non-zero value in the ordered phase. The order parameter is not unique, as its 
definition depends on the particular problem. In the Ehrenfest approach, the order of the transition is related to the 
order of the first discontinuous derivative of the energy. 
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FIG. 4. Order parameter (upper panel) and 2 + energy (lower panel) as a function of x for a system of N = 1000 bosons in 
the 0(6)-f/(5) transition. 



The advantage of an exactly solvable model is that we can reconcile these two approaches to phase transitions, 
by defining the order parameter as the derivative of the hamiltonian (6) with respect to x. Making use of the 
Helmann-Fcyman theorem, the order parameter is given by 

0= d -^=n d -±{S.d>-sU^){d.d-ss) - (0) = g 

In Fig. 4 we show the order parameter (we have subtracted the constant N — 1 to assure it is zero in the disordered 
phase) and the energy of the first excited state with angular momentum 2 + , for a system of N = 1000 bosons and in 
the vicinity of the phase transition. From fig. 4, we see that the phase transition is of second order. By including finite 
size corrections, we can readily study the detailed properties of this phase transition, as well as the corresponding 
phase transitions in larger spaces. The present analysis based on the exact solution for very large systems confirms 
recent results [11] based on the phenomenological Landau theory on the characteristics of this phase transition. 
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